Scattering lengths of calcium and barium isotopes 
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We have calculated the s-wave scattering length of all the even isotopes of calcium (Ca) and 
barium (Ba), in order to investigate the prospect of Bose- Einstein condensation (BEC). For Ca we 
have used an accurate molecular potential based on detailed spectroscopic data. Our calculations 
show that Ca does not provide other isotopes alternative to the recently Bose condensed 40 Ca that 
suffers strong losses because of a very large scattering length. For Ba we show by using a model 
potential that the even isotopes cover a broad range of scattering lengths, opening the possibility of 
BEC for at least one of the isotopes. 

PACS numbers: 34.20.-b, 34.50.Cx, 67.85.-d 



I. INTRODUCTION 

Knowledge of the s-wave scattering length plays a cru- 
cial role in the achievement of Bose-Einstein condensa- 
tion (BEC) in ultracold atomic gases A positive sign 
of the scattering length indicates that a Bose-Einstein 
condensate is stable whereas a negative sign indicates an 
unstable or collapsing condensate. Also its magnitude 
is crucial in the formation of BEC, as it determines the 
clastic and inelastic collision rates. For efficient evapora- 
tive cooling a large scattering length is required, however 
for a too large scattering length three-body recombina- 
tion loss limits the formation and collisional stability of 
a BEC. This gives rise to a range of about 50 — 200ao 
(ao=0.5292 nm) that can be considered favorable. 

Alkaline-earth elements and the alike ytterbium (Yb) 
system have been the center of attraction in ultra-cold 
atom research because of their unique atomic struc- 
ture and the availability of different stable isotopes with 
zero nuclear spin (bosons) and with large nuclear spin 
(fermions). They are being used in frequency metrol- 
ogy @~S I an0 - are proposed for quantum information pro- 
cessing 5|. Recently, BEC has been achieved in 40 Ca @, 
84 Sr @,H, 86 Sr [| and 88 Sr [Tjj. Also BEC of different 
Yb isotopes has been reported: 174 Yb [ill, 170 Yb [Hj], 
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Yb [13J and Yb [14[, of which the last one has an 
abundance of only 0.13%, highlighting the possibility of 
using rare isotopes to achieve BEC. 

In general, isotopes of an element have different scat- 
tering lengths. Therefore the ability of BEC formation 
depends crucially on the chosen isotope. For instance, 
the scattering length of 84 Sr is 123ao [la]; which made it 
an ideal candidate to achieve BEC despite its low abun- 
dance of only 0.6% fl6| . as compared to the more abun- 



dant 86 Sr (10%, 800a [H) and 88 Sr (83%, -2a [Hj), 
either suffering strong losses |9| or requiring sympathetic 
cooling with another isotope respectively. The 40 Ca 
isotope has a large scattering length of about 440ao [|[, 
limiting the size and stability of BEC. 

Motivated by the success of BEC in Sr and Yb even 
isotopes, and 40 Ca, we carried out calculations to obtain 
the scattering lengths for all the even (bosonic) isotopes 
of Ca and Ba, for which mass and abundance are summa- 
rized in Table |U Laser cooling and trap ping of all stable 
Ca isotopes have been reported [ij], [Hf ■ The demonstra- 
tion of magneto-optical trapping of Ba [lj| has opened 
its use for ultra-cold collision studies, photoassociation 
spectroscopy and Bose-Einstein condensation. Recently, 
an optical clock based on ultracold Ba has been pro- 
posed 0. 



II. THEORY 

The scattering properties of ground state atoms are 
obtained from the underlying two-body ground state po- 
tentials. For alkaline-earth atoms there is only one molec- 
ular ground state potential, namely the singlet A^E^" po- 
tential. Furthermore, the lack nuclear spin and therefore 
hyperfine structure for the even (bosonic) isotopes highly 



TABLE I: Mass and abundance of all stable even isotopes of 
Ca and Ba El. 
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Calcium 


Isotope 


40 


42 


44 


46 


48 


Mass (amu) 


39.963 


41.959 


43.955 


45.954 


47.953 


Abundance (%) 


96.941 


0.647 


2.086 


0.004 


0.187 


Barium 


Isotope 


130 


132 


134 


136 


138 


Mass (amu) 


129.906 


131.905 


133.904 


135.905 


137.905 


Abundance (%) 


0.11 


0.10 


2.42 


7.85 


71.70 
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FIG. 1: (Color online) The X X T,+ potentials of Ca 2 and Ba 2 
as a function of internuclear distance. The solid blue curve 
is the potential of Ca 2 taken from Ref. [53]. The dashed 
red curve is the Tang-Toennies potential of Ba 2 taken from 
Ref. Hi. 



FIG. 2: (Color online) The mass dependence of the scattering 
length for Ca, using the X 1 ^ potential of Ref. [23|], showing 
calculations with Ca coefficients that give rise to 340ao (red 
closed circles) and 700ao (blue closed squares) for 40 Ca. The 
plotted mass is twice the reduced mass and the dashed vertical 
lines indicate the mass of the even isotopes of Ca. 



reduces the number of collision channels compared with 
the alkali metal systems. 

To calculate the s-wave scattering length we solve the 
ID radial Schrodinger equation with zero angular mo- 
mentum and vanishing kinetic energy, 



2f "\F 



V(R)]i>(R) = Q, 



(1) 



where \i is the reduced mass, R the internuclear dis- 
tance, V(R) the Born-Oppcnheimer potential (here the 
X : E+ potential), and E the kinetic energy (below fiK). 
The asymptotic form of the wavefunction is ip(R) oc 
sin[fc(i? — a)], where k = ^/2^iE/h 2 , and a is the scatter- 
ing length. We fit the asymptotic form to the solution of 
Eq. (JTJ for large R in order to obtain a. Within the Born- 
Oppcnhcimcr approximation, the potentials arc identical 
for all the isotopes of a particular atomic system. There- 
fore, for a given potential V(R), one only needs to adjust 
/i to obtain a for all the isotopes. 

To obtain reliable scattering lengths accurate poten- 
tials are needed. With the exception of systems with 
only a few electrons, like metastable helium (22[, ab ini- 
tio potentials are in general not accurate enough and con- 
straints from experimental data are required. For Ca 2 an 
analytic representation of the potential based on 

an extensive set of experimental data from the Tiemann 
group is available [23| . For Ba2 such an accurate po- 
tential is not available. Here we rely on the analytical 
representation according to the Tang-Toennies potential 
model 24 1, which has shown to be able to reproduce ac- 



III. RESULTS 



Calcium 



curately the X X T,+ potentials of Ca 2 [25| and Sr 2 [2^ 
The used Ca 2 and Ba 2 potentials are shown in Fig. [1] 



Information on the Ca 2 potential has been 

gathered by several spectroscopic methods, including 
photoassociation spectroscopy [26ti28l | , Fourier-transform 
spectroscopy [2!| and filtered laser excitation tech- 
nique [23[. Based on these detailed spectroscopic data, 
an interval of the scattering length for 40 Ca is determined 
to be 340 - 700a [H|. In the 40 Ca BEC experiment a 
scattering length of a « 440eto was estimated from a mea- 
surement of the chemical potential 0, without quoting 
an uncertainty. Therefore we take the scattering length 
range of 340 — 700ao as starting point of our calculations. 

We have taken the analytical representation and pa- 
rameter values of the Ca 2 X 1 Y<+ potential from Ref. [23| 
(see Fig. [T]). We allow the Cq long-range coefficient to 
vary within a factor 0.991 to 1.003 from its reported value 
[23[, in order to reproduce the scattering length interval 
for 40 Ca. We then take the mass as a variable parame- 
ter and calculate the scattering length for the full mass 
interval spanned by the stable even isotopes. In this way 
we transfer the knowledge on 40 Ca to the other isotopes. 

The results for Ca are shown in Fig. [21 The plotted 
mass is twice the reduced mass, which for homonuclear 
collisions simply is the atomic mass. The dashed vertical 
lines indicate the masses of the even isotopes, from which 
the scattering lengths for homonuclear collisions are di- 
rectly read off. Two calculations are shown, obtained 
from potentials with C§ coefficients that give rise to 340ao 
(red closed circles) and 700ao (blue closed squares) for 
40 Ca. The scattering length shows the expected behavior 
as a function of mass, with a regular pattern of scatter- 
ing resonances. They can be understood from the mass 
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TABLE II: Scattering lengths of even Ca isotopes (left), as 
well as interisotopic scattering lengths for all combinations 
with 40 Ca, and 42 Ca+ 44 Ca (right). 



a (ao) 


a (ao) 


40+40 


+340 . . . + 700 [28] 


40+42 


+45 . . . + 50 


42+42 


+320 . . . + 640 


40+44 


±oo 


44+44 


+460 . . . + 1800 


40+46 


+62 . . . + 68 


46+46 


±oo 


40+48 


-66 39 


48+48 


-230... -120 


42+44 


+47 . . . + 52 



dependence of the vibrational splitting, and therefore the 
number of bound states. A scattering resonance appears 
at those mass values at which a new vibrational state 
becomes bound, i. e. where the least bound vibrational 
state has zero binding energy. We find that accidentally 
all even isotopes are located close to such a scattering 
resonance, giving rise to large positive or negative scat- 
tering lengths. 

The scattering length intervals for the different iso- 
topes are given in Table |TIJ It is clear that Ca does not 
provide an isotope with a favorable scattering length. At 
most one can state that 42 Ca has a slightly smaller scat- 
tering length than 40 Ca. In addition, the intcrisotopc 
scattering lengths are given for collisions between 40 Ca 
and the other isotopes, and 42 Ca+ 44 Ca. Here we find 
favorable scattering lengths for 40 Ca+ 42 Ca, 40 Ca+ 46 Ca 
and 42 Ca+ 44 Ca. Taking into consideration the large scat- 
tering length of 44 Ca and 46 Ca, the only interesting mix- 
ture in view of sympathetic cooling is 40 Ca+ 42 Ca. 



B. Barium 

Knowledge of the Ba2 X 1 ^ potential is sparse, be- 
cause of a limited amount of experimental data |30L l3lT ] , 
which cover only a small part of the vibrational spectrum. 
Similarly, theoretical data is scarce. The binding energies 
of B&2 were calculated by Jones [HI]. First wave func- 
tion based quantum chemical calculation was reported 
by Ref. [HI, most recently by Ref. [34|. The most ac- 
curate calculations of the long-range coefficients Cr, Cg, 
Cio are reported by Porsev and Derevianko [35l l36| . 

We have used the Tang-Toennies potential [331 of 
Ref. [Hi to model the Ba 2 X 1 ^ potential as 



V(R) = Ae 



-bR 



E 



-bR 



E 

k=0 



(buy 




(2) 



where the short-range parameters (in a.u.) A = 105.4 
and b = 0.9657 are obtained from rescaling an accurate 
Si"2 Ar x S+ potential (38l |. using the measured fundamen- 



tal vibration frequenc 
coefficients from Ref. 135 



3 If. and taking the long-range 
~36| (see Fig.p. 
Although we do not expect the potential to be accurate 
enough to predict the scattering length, we can use it to 
investigate its mass dependence. The result is shown in 
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FIG. 3: (Color online) The mass dependence of the scattering 
length for Ba, using the Tang-Toennies model to represent 
the -X^El" potential with the parameter values from Ref. [24[ 
(black closed circles). The arrow indicates the effect of the 
uncertainty in the theoretical Ca coefficient, showing its lower 
(blue closed triangles) and upper (red closed squares) limit. 
The plotted mass is twice the reduced mass and the dashed 
vertical lines indicate the mass of the even isotopes of Ba. 



Fig. |3] (black closed circles). We observe that the spac- 
ing between the resonances is larger than that between 
the even isotopes, giving rise to a variety of scattering 
lengths among the different isotopes. Therefore it is very 
probable that at least one of the isotopes has a favorable 
scattering length for BEC formation. 

In addition, we have investigated the effect of the 
theoretical uncertainty in the reported Cg coefficient, 
5160 ± 74 a. u. [35ll3^ |. We have repeated the above cal- 
culation with the upper (red closed squares) and lower 
(blue closed triangles) limits (see Fig. [3]). Note that 
these a djust ments in Cq also require small changes in 
A and b [2J]. The uncertainty in Cq gives rise to a large 
spread in the scattering length, which can also be seen 
as a large shift along the mass axis, as indicated by the 
arrow. This means that even if an accurate short-range 
potential would be available, still no precise prediction 
of the scattering length is possible. On the other hand, 
whereas the scattering length itself is extremely sensi- 
tive to details of the potential, the spacing between the 
scattering resonances is not. Therefore, once scattering 
length information of one of the isotopes becomes avail- 
able, the simple Tang-Toennies potential is sufficient to 
predict the scattering lengths of all the other isotopes 
[40| . at least at an accuracy that allows to choose the 
best isotope to achieve BEC. 



IV. CONCLUSION AND OUTLOOK 

In summary, we have investigated the s-wave scatter- 
ing lengths of the Ca and Ba even (bosonic) isotopes. An 
experimentally accurate ground state potential of Ca 2 
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enabled us to calculate the scattering lengths for all the 
isotopes, taking knowledge of 40 Ca scattering length as a 
reference. We have observed that accidentally scattering 
resonances appear near all even isotopes, leading to large 
positive and negative scattering lengths for all isotopes. 
Only 42 Ca shows a small improvement in terms of scat- 
tering length compared to 40 Ca. Therefore, alternative 
methods to that of evaporative cooling arc interesting for 
Ca, such as direct laser cooling into quantum degeneracy 
using a narrow transition (39j . 

In addition, scattering lengths for all the even Ba 
isotopes have been investigated using a Tang-Toennies 
model potential. In contrast to the Ca case, here the 
scattering lengths vary strongly over the isotopes, result- 
ing in a high probability that at least one isotope has 
a favorable scattering length to produce a BEC. More 
accurate knowledge of the ground state Ba2 potential, 
including the long-range coefficients, based on experi- 
mental data that covers the full vibrational spectrum is 



needed for an improved determination of the scattering 
lengths. However, once experimental knowledge on the 
scattering length of one isotope becomes available, the 
scattering lengths of all the other isotopes can be pre- 
dicted with the present model potential. 
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